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0> ' We study the effects of the Coulomb interactions between electrons on the 

Hofstadter butterfly, which characterizes the subband structure of the Landau 
levels of a two-dimensional electron gas in a perpendicular homogeneous magnetic 
<^ ! field and a periodic lateral superlattice potential. The interactions essentially 

preserve the intricate gap structure of the Hofstadter spectra, but with a lower 
symmetry that depends on the filling of the Landau bands. For short enough 
periods and strong enough modulation the miniband structure can be resolved 
^ . in the thermodynamic density of states. 

^ ■ Measurements of the magnetotransport in superlattices weakly modulated in one or two 

OO ' directions have revealed novel oscillations in the conductivity, the Weiss oscillations, reflect- 

^ . ing the commensurability of the cyclotron radius Rc = I'^kp = vp/uj^ of the electrons at 

Q> ! the Fermi energy Ep = fi^kp/ {2m*) and the modulation period L; I = {ch/eBY^'^ is the 

magnetic length, and Uc = eB/{mc) the cyclotron frequency IQ. The Weiss oscillations 
g ' are superimposed on the Shubnikov-de Haas oscillations which are well known in a homo- 

geneous two-dimensional electron gas (2DEG) and caused by the commensurability of the 
Fermi wavelength \p and /. The single-particle energy spectrum for this problem has been 



O . investigated by several researchers culminating in Hofstadter's butterfly, a graph showing 

^ ' the complicated self-similar splitting of a Landau band (LB) into minibands as a function 

of the magnetic flux through a unit cell of the lattice (for a review see |@]). The effects of 
^ ■ the splitting of LB's into minibands has, to the best of our knowledge, not been observed 

■ directly in the currently common superlattices (L ~ 300 nm = 3 x 10^^ cm~^) since usu- 

ally several LB's are occupied, which may overlap near Ep. The energy spectrum has been 
investigated in the intermediate region in the absence of collision broadening. The coupling 
of the LL's by the external periodic potential strongly reduces the original high symmetry 
of the Hofstadter spectrum but retains a very comphcated subband structure 0, ^. All the 
above-mentioned theoretical investigations of the energy spectrum of a 2DEG in a periodic 
potential and a homogeneous perpendicular magnetic field have neglected the effects of the 
electron-electron interaction. Here we study, within the Hartree approximation (HA), these 
effects on the electronic energy spectrum and the thermodynamic density of states. 

To describe the electrons in the conduction band in the presence of a lateral superlattice 
at the AlGaAs-GaAs interface in a constant perpendicular magnetic field B = Bz we employ 
a model of a 2DEG, with the three-dimensional charge density given by —ens{r)5{z), and 
r = {x,y). The square superlattice is spanned by orthogonal lattice vectors of length L, 
the primitive translations of the Bravais lattice B. The corresponding reciprocal lattice TZ 
is spanned by G = Gigi + (^2^2, with Gi,G2 E Z, gi = 2nx/L, and ^2 = 27ry/L. The 
external periodic potential the electrons are moving in is taken to be of the simple form 
V{r) = Vo{cos(5'ia;) -|- cos{g2y)}. The electron-electron interaction is included in the HA 
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leading to an effective single-electron Hamiltonian H = Hq + Vnir) + V{'r), where Vh(^) is 
the effective potential in a medium with a dielectric constant k, felt by each electron and 
caused by the total charge density of the 2DEG, —ens{r), and the neutralizing background 
charge density +enb = +e{ns{r)). The periodic external potential V{r) and the constant 
external magnetic field imply that all physical quantities of the noninteracting system are 
periodic with respect to translations oi R E B. The periodicity of the Hartree potential 
follows from that of Usif). In the following calculations we use the symmetric basis functions 
constructed by Ferrari and used by Silberbauer . The set of Hartree equations has to be 
solved iteratively together with the condition that the average electron density = Ng/A 
is constant, which determines the chemical potential /i. Ng is the number of electrons per 
unit cell with area A in B. 

The Coulomb interaction not only couples directly the subbands of a particular LL, but 
also the subbands of different LL's. If only the basis states in the same LL are taken into 
account, the modulation strength is irrelevant since it can be factored out of the Hamiltonian 
matrix. In this situation, corresponding to the usual discussion of Hofstadter's butterfly 
on the basis of Harper's equation, a highly symmetric energy spectrum is expected. Of 
course, this restricted model is only appropriate for describing a system with very weak 
modulation (as compared with the inter-Landau-level energy spacing huc)- The bandwidth 
of the subbands of the first LB is shown in Fig. |l], for the two cases of 3 or 6 magnetic 
flux quanta through the unit lattice cell, •I'/'I'o = P"? = 3, 6, identifiable from the number of 
subbands appearing in each figure. The period is short, L = 50 nm, the modulation strong, 
Vq = 4.0 meV, and T = 1.0 K. However only one LL is used in the calculation, so that the 
results apply to the case of a weak superlattice potential. The electron-hole symmetry is 
exact, as is reflected in each case by the point symmetry of the subband structure around 
fi at half filling. The screening is strongly dependent on Ng and the average filling factor. 
It is, therefore, strongest for pq = 6 when maximum six electrons can occupy the lowest 
Landau level in the unit cell and when the LB is half filled, A^^ = 3. The (dimensionlessly 
written) thermodynamic density of states Dt = t^TiuJc{dns/ dix)TBi with Hs = {ns{r)), is 
shown in Fig. |[ Clearly the strong screening in the pq = 6 case obscures the fine structure 
in Dt leading to only one pronounced peak. The subbands in the other cases turn up as 
separate peaks in Dt- The Hofstadter butterfly should be recovered when the calculation is 
performed for A^^ = and only one LL, and if the width and location of the subbands are 
plotted as functions of the inverse of the number of magnetic flux quanta through a unit cell 
of the lattice, i.e. of $o/^ = addition, all energies have to be scaled according 

to: E ^ {E — hujc/2) /Vq. For interacting electrons A^^, the average filling factor of the LB 
turns out to be a new parameter controlling the bandstructure. We therefore, present in Fig. 
1^ the subband structure for four values of the inverse flux l/{pq) in the case of A^^ = 1-50 
together with the location of fi. For low density, the subband structure is almost symmetric 
around the energy zero, like in the noninteracting case, but the subbands become quite 
asymmetric for a higher density of electrons. For the short-period superlattice studied here, 
and for a low density of electrons, the essential gap structure does survive in the presence of 
interaction. The presentation in Fig. ^ corresponds to the experimental procedure of keeping 
the density of electrons fixed but changing the magnetic field. Another way to investigate 
the screening is to keep the average filling factor u constant but vary the magnetic field 
and, thus, also the average density of the 2DEG. Fig. |^ compares the subband structures for 
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= l/2j 1/3) 1/4) 1/6 and v = 1/2 with the complete Hofstadter butterfly. Here the 
energies have been scaled with the factor {E — huJc/2) exp{(7rZ/L)^}/Vo so that the results can 
be directly compared with earlier calculations of the Hofstadter spectrum [Q]. The energy 
spectrum in the interacting case shows an overall reduction in dispersion or bandwidth due 
to the strong screening that is most effective for large flux and large number of available 
states. The bandwidths of the subbands for the interacting case has been evaluated here on a 
discrete lattice in the magnetic Brillouin zone without attempting an interpolation between 
the lattice points; thus, the actual bandwidths can be larger by a small percentage of the 
widths shown. 

In an interacting 2DEG subject to a superlattice potential and a homogeneous perpen- 
dicular magnetic field not only the magnetic flux through a unit cell but also the density of 
electrons determines the complicated splitting of the LL's into subbands. We have shown 
that in the HA the essential gap structure of the energy spectrum remains, although the 
screening leads to a quenching of the Hofstadter butterfly at small values of the inverse flux. 
The symmetry of the butterfly is also lowered by the coupling to higher LL's due to increased 
strength of the periodic potential, as has already been discussed by other authors ^. For 
periods around L = 200 nm (currently attainable in experiments on superlattices) the 2DEG 
can effectively screen the periodic potential even for a very low density Ug when only one 
or two Landau bands are partly occupied. Only at shorter lattice constants (L < 100 nm) 
and thus much higher magnetic fields we can, on the basis of our Hartree calculation, expect 
the subband structure to be resolvable in experiments, when Ug is maintained low enough. 
These predictions are made on the basis of the calculated energy spectrum and the struc- 
ture observed in the thermodynamic density of states Dt- The transport coefficients do not 
depend on in any simple way, and may be more sensitive to the subband structure of the 
energy spectrum . In order to keep the computational efforts in reasonable limits, we have 
restricted our calculations to a few characteristic values of the magnetic field. Our results in- 
dicate in which manner screening effects will change the overall appearence of the Hofstadter 
energy spectrum. Unfortunately, there is no simple way to extrapolate from these special 
values to arbitrary rational numbers of flux quanta per unit cell. From the experimental 
point of view these details of the energy spectrum seem not to be accessible anyway. The 
challenge for the near future is to resolve the most prominent gaps of the energy spectrum 
experimentally. 
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Figure 1: The bandwidth of the minibands of the lowest Landau band as a function of Ng for 
a Hartree interacting 2DEG and {p — 1, q — 3) (left), (p = 2, g = 3) (right). The chemical 
potential fi is indicated by a continuous curve. L = 50 nm, T = 1.0 K, Vq = 4.0 meV. One 
Landau level is included in the HA. m* = 0.067mo, k = 12.4. 



Figure 2: The thermodynamic density of states (left) Dt = PhuJc{dns/d^)TB as a function 
of 1/pq and Ns. The scaled bandwidth (right) of the subbands of the lowest Landau band, 
{E — ^a;c/2)/Vo, as function of l/pq for Ns = 1.5. The chemical potential is indicated by 
crosses. L = 50 nm, T — 1.0 K, Vq ~ 4.0 meV. One Landau level is included in the HA. 



Figure 3: The scaled bandwidth {E — huj,,/2) exp{(7r//L)^}/Vo of the subbands of the lowest 
Landau band as function of l/pq for the Hartree interacting 2DEG with u = 1/2 (left), and 
for the noninteracting 2DEG (right), the Hofstadter butterfly. One Landau level is included 
in the calculations. In the left subfigure the 6 minibands for the case pq = 6 can not all 
be resolved due to vanishing band gaps. The parameters used for the left subfigure are: 
L = 50 nm, T = 1.0 K, Vq = 4.0 meV. 
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